Low energy spin waves and magnetic interactions in SrFe 2 As 2 
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We report inelastic neutron scattering studies of magnetic excitations in antiferromagnetically 
ordered SrFe2As2 (Tjv = 220 K), the parent compound of the FeAs-based superconductors. At low 
temperatures (T = 7 K), the spectrum of magnetic excitations S(Q,hui) consists of a Bragg peak 
at the elastic position (hu = meV), a spin gap (A < 6.5 meV), and sharp spin wave excitations 
at higher energies. Based on the observed dispersion relation, we estimate the effective magnetic 
exchange coupling using a Heisenberg model. On warming across Tn, the low temperature spin gap 
rapidly closes, with weak critical scattering and spin-spin correlations in the paramagnetic state. 
The antiferromagnetic order in SrFe2As2 is therefore consistent with a first order phase transition, 
similar to the structural lattice distortion. 



The parent compounds of the high-transition tempera- 
ture (high-T c ) copper oxides are simple antiferromagnetic 
(AF) Mott insulators [H characterized by a very strong 
nearest neighbor AF exchange coupling J (> 100 meV) in 
the CuC>2 planes Q . When holes or electrons are doped 
into the Cu02 planes, the character of the ground state 
is fundamentally altered from an AF insulator to a su- 
perconductor with persistent short-range AF spin corre- 
lations (excitations) Q . In the case of FeAs-based super- 
conductors such as iJFeAsOi-^Fa; (where R = La, Nd, 
Sm, Pr,...) 0, H| i, and A 1 _ x B x Fe 2 As 2 (A =Ba, Sr, 
Ca, B = K, Cs, Na) H, H E3, E2 , although the undoped 
parent compounds arc also long-range ordered antifer- 
romagnets with a collinear spin structure as shown in 
Fig. la [H O, 0, [H, EE E3, El, much is unknown 
about the magnetic exchange coupling responsible for 
such a spin structure. For example, early theoretical 
studies suggested that LaFeAsO has a spin-density-wave 
(SDW) instability [l§ El- As a consequence, the AF 
spin structure in these materials arises from quasiparti- 
cle excitations across electron-hole pockets in a nested 
Fermi surface 2l|, much like SDW antiferromagnetism 
in metallic chromium (Cr) (22j. Alternatively, a Heisen- 
berg magnetic exchange model[23l. liil. |25| is suggested to 
explain the AF structure. Here, the collinear spin phase 
is stable when the nearest neighbor exchange J\ and the 
next nearest neighbor exchange J 2 satisfy J\ < 2 J 2 (Fig, 
la). First-principles calculations estimate J\ ~ Jo j24| . 
In contrast, some band structure calculations [26[ suggest 
that the J\ along the a-axis and &-axis of the low tem- 
perature orthorhombic structure (c > a > b) can have 
different signs with J\ a and Ju being AF and ferromag- 
netic, respectively, and that J\ a > 2J 2 . Therefore, there 
is no theoretical consensus on the relative strengths of 
Jia, Ju, and J 2 or the microscopic origin of the observed 
AF spin structure. If magnetism is important for super- 
conductivity of these materials, it is essential to establish 



the "effective Hamiltonian" that can determine the mag- 
netic exchange coupling. 

In this paper, we report inelastic neutron scattering 
studies of spin wave excitations in single crystals of 
SrFc 2 As 2 (T N = 220 K) [T^]. At low temperature, we 
find that spin waves have an anisotropy gap of A = 6.5 
meV and disperse rapidly along both the [H, 0, 0] and 
[0, 0, L] directions. On warming to 160 K, the magnitude 
of the spin gap decreases to 3.5 meV while the intensi- 
ties of the spin wave excitations follow the expected Bosc 
statistics. However, there are only weak critical scatter- 
ing and magnetic correlations in the paramagnetic state 
at 240 K, in sharp contrast to the SDW excitations in 
Cr [2H and spin waves in cuprates 0, Q ■ We estimate 
the effective magnetic exchange coupling using a Heisen- 
berg model and find that J\ a + 2J 2 = 100 ± 20 meV, 
J z = 5 ± 1 meV, with a magnetic single ion anisotropy 
J s = 0.015 ±0.005 meV. The weak critical scattering and 
paramagnetic spin-spin correlations, together with the si- 
multaneous first order structural transition p~7i | . suggest 
a first order AF phase transition. 

Our experiments were carried out on HB-1 triple- axis 
spectrometer at the High Flux Isotope Reactor, Oak 
Ridge National Laboratory, BT-7 and SPINS triple-axis 
spectrometers at the NIST Center for Neutron Research. 
For HB-1 and BT-7 measurements, we fixed the final 
neutron energy at Ef = 14.7 meV and used PG(0,0,2) 
(pyrolytic graphite) as monochromator and analyzer. A 
PC filter was placed in the exit beam path to eliminate 
A/2. For the SPINS measurements, the final neutron en- 
ergy was fixed at Et = 4.9 meV and a cold Be filter 
was placed in the scattered beam path. SrFc2As2 single 
crystals were grown from flux Q and coaligned within 
2 degrees to have a total mass of ~0.7 g. From earlier 
diffraction work 17|, we know that the AF order occurs 



in close proximity to the lattice distortion, changing the 
crystal structural symmetry from tetragonal above TV 
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FIG. 1: (Color online) (a) The Fe spin ordering in the 
SrFe2As2 chemical unit cell and magnetic exchange couplings 
along different high-symmetry directions, (b) The AF Neel 
temperature and the temperature dependence of the struc- 
tural (2, 2, 0) Bragg peak of the SrFe 2Sr2 crystals used in the 
experiment [lj. The inset shows positions in reciprocal space 
probed in the experiment, (c) Observed spin wave dispersion 
along the [H, 0, 0] direction at 160 K. (d) Similar dispersion 
along the [0, 0, L\ direction, (e) Calculated three-dimensional 
spin wave dispersions using Ji a = 20, Jib = 10, Ji = 40, 
J z — 5, and J s — 0.015 meV. (f) Temperature dependence of 
the anisotropy spin gap A(T). 



to orthorhombic below it (Fig. lb). However, it is un- 
clear whether the structural and ma gnet ic phase tran- 
sitions are second [27| or first order [28j . For the ob- 
served spin structure (Fig. la), magnetic Bragg peaks 
are allowed at [H, 0, L] (H, L are odd integers) recipro- 
cal lattice units (r.l.u), where the momentum transfer is 
Q(in A" 1 ) = (H2n/a, K2ir/b, L2n/c) and a = 5.5695(9), 
b = 5.512(1), c = 12.298(1) A are lattice parameters in 
the orthorhombic state at 150 K. To probe spin wave exci- 
tations, we aligned our single crystal array in the [H, 0, L] 
zone, where we can probe excitations along the [H, 0, 0] 
and [0,0, L] directions. 

Figures 2a and 2e show constant-energy scans for 
hu) = 1 and 5 meV around [H, 0, 1] at 160 K obtained on 
SPINS. While the scattering at hcu = 1 meV is feature- 
less (Fig. 2a), there is a clear peak centered at H = 1 
in the 5 meV data (Fig. 2e). This immediately sug- 
gests that spin waves in SrFe2As2 have an anisotropy 
gap at this temperature that is less than 5 meV. Mov- 
ing on to higher energies, Figures 2b-d and 2f-h summa- 




FIG. 2: Wave vector dependence of the spin wave excita- 
tions at 160 K obtained on cold [(a) and (e)] and thermal 
[(b)- (d) and (f)-(h)] triple-axis spectrometers at different en- 
ergies, (a) Q-scan along the [H, 0, 1] direction at fwj = 1 meV 
using SPINS. The spectrum is featureless indicating the pres- 
ence of a spin gap exceeding 1 meV. Identical scan at Tilu — 5 
meV in (e) shows clear evidence of spin wave excitations cen- 
tered at (1, 0, 1). (b-d) Q-scans along the [H, 0, 1] or [H, 0, 3] 
directions at different energies. The spectra clearly broaden 
with increasing energy, (f-h) Similar scans along the [1,0, L] 
direction, which probe the exchange coupling J z . 



rize Q-scans along [H, 0, 0] and [0, 0, L] directions, respec- 
tively, at different energies. The Q-widths of the scatter- 
ing clearly become broader with increasing energy. Fig- 
ures lc and Id show the observed dispersion curves for 
the limited energy range with observable spin wave exci- 
tations. Assuming an effective Hciscnbcrg Hamiltonian 

H El h = j la e 4 s t ■ s s + J u E< Si • Sj + -h T,a s, • 



Sj + J z Si ■ Sj - J s (Sf) 2 , where J la , J u , J 2 , and J z are 
exchange interactions shown in Fig. la; J s is the sin- 
gle ion anisotropy; and S is the magnitude of iron spin, 
the spin wave dispersions along the [H, 0, 0] and [0, 0, L] 
directions near the (1,0,1) Bragg peak are E(k x ) = 
2S[{J la + 2J 2 + J S + J z f - (J z - (J lQ + 2J 2 ) cosK) 2 ] 1 ' 2 
and E(k z ) = 2S\(2Ji a + 4 J 2 + J s +J Z -J Z cos k z )(J s + 
J z + J z cos kz)] 1 / 2 , respectively. In addition, the size 
of the spin gap due to the single ion anisotropy is 
A(l, 0, 1) = 2S[J s {2J la + 4J 2 + J s + 2J Z )\ X ' 2 . The solid 
lines in Figs, lc and Id are the best fits with these equa- 
tions, where J la + 2J 2 = 100 ± 20 meV, J z = 5 ± 1 meV, 
and J s = 0.015 ±0.005 meV. The three-dimensional plot 
in Fig. le shows the expected spin wave dispersion at 
higher energies. 

To determine the temperature dependence of the spin 
gap, we carried out energy scans at the signal and back- 
ground positions for spin waves at different temperatures. 
At 7 K, an energy scan at the magnetic zone center po- 



3 




10 12 14 16 



FIG. 3: (Color online) Temperature dependence of the spin- 
gap obtained from energy scans around the (1,0,1) and (1,0,3) 
Bragg peaks, (a) Low temperature (T = 7 K) constant- 
ly scans at the signal (Q — (1,0,1)) and background (Q = 
(1.2,0,1)) positions show a clear spin-gap of A = 6.5 meV 
(b) Similar scans at Q = (1,0,3) and Q = (0.8,0,3) which 
again show A = 6.5 meV. (c),(d) Temperature dependence of 
the spin-gap, where A = 4.5 meV at 80 K and A = 3.5 meV 
at 160 K. 



sition Q = (1,0,1) shows an abrupt increase above 6.5 
meV, while the background scattering at Q = (1.2,0,1) 
is featureless (Fig. 3a). Energy scans at equivalent posi- 
tions Q = (1,0,3) and (0.8,0,3) in Fig. 3b show similar 
results and therefore reveal a low temperature spin gap 
of A = 6.5 meV. On warming to 80 K, identical scans at 
Q = (1,0, 1) and Q = (1.2, 0, 1) show that the spin gap 
is now at A = 4.5 meV (Fig. 3c). Finally, A becomes 
3.5 meV at 160 K, consistent with constant-energy scans 
in Figs. 2a and 2e. These results indicate that the spin 
anisotropy of the system reduces with increasing temper- 
ature. 

Figure 4 summarizes the temperature dependence of 
the spin waves and quasielastic scattering in the AF or- 
dered and paramagnetic states obtained on SPINS and 
BT-7. Figure 4a shows energy scans at Q = (1, 0, 1) for 
T = 130, 160, 200, 240, and 280 K plotted on a log scale. 
When the temperature is increased across Tn (= 220 K), 
there is a rapid decrease in the ordered moment but lit- 
tle evidence for quasielastic and critical scattering, which 
are signatures of a second order phase transition. To il- 
lustrate this point, we plot in Fig. 4b the temperature 
difference scattering using 280 K data as background. 
Besides the magnetic Bragg peak below Tn at Tiuj = 
meV, there is little quasielastic critical scattering typi- 
cal of a second order phase transition. Figure 4c shows 
constant-energy scans (hu) = 1 meV) measured on SPINS 
and the scattering is essentially featureless at all temper- 
atures investigated. Assuming no spin correlations in the 
paramagnetic state at 280 K, the differences in the scat- 
tering between 240 K and 280 K at Tiui = 1 meV should 
reveal the magnetic intensity gain close to Tn- Con- 
sistent with the temperature dependence of the energy 
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FIG. 4: (Color online) Temperature dependence of the quasi- 
elastic magnetic scattering and spin wave excitations below 
and above Tn- Data in (a)-(c) are obtained on SPINS, (a) 
Constant-Q scans at the Q = (1,0, 1) Bragg peak position at 
different temperatures. Except for the dramatic increase in 
the elastic component below Tn, the quasi-elastic scattering 
above hui = 0.5 meV is essentially temperature independent, 
revealing no evidence for the Lorentzian-like paramagnetic 
scattering above Tn observed in Cr [2^]. (b) Temperature 
difference spectra using T = 280 K scattering as background. 
The data again show little evidence of critical scattering at 
Tiuj — 1 meV. (c) Q -scans at hui = 1 meV at different tem- 
peratures. We speculate that the slight increase in overall 
scattering at 240 K from 280 K shown in (e) is due to weakly 
correlated paramagnetic spins, (d) Tilu — 10 meV spin-wave 
excitations at 160 K and 200 K obtained on BT-7. (f) hui = 16 
meV spin-wave excitations at 7 K and 160 obtained on HB-1. 
The intensity increase is due to the Bose population factor. 



scans in Figs. 4a and 4b, there are signs of possible un- 
corrected paramagnetic scattering (since the subtracted 
data in Fig. 4e are overall positive) at 240 K but weak 
critical scattering. For temperatures below Tn, we find 
that spin wave excitations at temperatures below 160 K 
simply follow the Bose statistics (Fig. 4f). 

The discovery of the collinear AF order with small mo- 
ment in LaFeAsO 12| has caused much debate about 
its microscopic origin. Since LaFeAsO is a semimetal, 
the observed AF order may arise from a SDW instability 
due to Fermi surface nesting [l9|, [2(| 21 1, where elec- 



tron itinerancy is important much like incommensurate 
SDW order in pure metal Cr [22j ■ Alternatively, there 
are reasons to believe that LaFeAsO is in proximity to a 
Mott insulator [3(j, and the AF order is a signature of 
local physics and electron correlations 25j, |3l| . Another 



heavily debated issue is the first [28| or second [27J or- 
der nature of the simultaneous structural/magnetic phase 
transition in SrFc2As2- 
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If the observed AF order in SrFe2As2 originates from 
Fermi surface nesting similar to the SDW order in 
Cr, the velocity of the spin waves c should be c = 
y/ v e Vh/3, where v e and Vh are the electron and hole 
Fermi velocities, respectively 22|. The dispersion re- 
lation is then hu> = cq where q is the magnitude 
of the momentum transfer away from the Bragg posi- 
tion. For Cr, the spin-wave velocity is measured to 
be c = 851 ± 98 meV-AQ. In addition, there are 
strong spin-spin correlations in the paramagnetic state 
where the dynamic structure factor S(q, Hlo) can be de- 
scribed by the product of a Gaussian centered at the 
SDW ordering wave vector and a Lorcntzian in en- 
ergy, or S(q,hu) = So{T)e-^ 2 {huj/k B T)/[{{huj) 2 + 
r)(l — g-ftw/fcsT^ wnere a and T are the Gaussian and 
Lorentzian widths, respectively At temperatures 

as high as 500 K (T = 1.6TV), one can observe a clear 
resolution-broadened Lorentzian centered at hu = meV 
with r = 15.6 meV (Hf. For comparison, there is no 
evidence of a Lorentzian-like quasielastic scattering in 
SrFe2As2 even at T = 1.09XV. The lack of critical scat- 
tering both below and above T/v, together with the fact 
that there is also an abrupt structural distortion occur- 
ring at the same temperature and observed thermal hys- 
teresis 



17 . Til 27 . liil ]. is consistent with the AF phase 
transition being first order in nature. 

To compare the observed exchange couplings in Fig. 
1 and those expected from SDW excitations in a nested 
Fermi surface, we note that Fermi velocities estimated 
from the local density approximation calculations for 
BaFe 2 As 2 are v e = 2.2 eV-A and v h = 1.2 eV-A. 



Assuming BaFc2As2 and SrFe2As2 have similar Fermi 
velocities, the expected spin wave velocity is then c ~ 
0.94 eV-A. However, since Angle Resolved Photoemis- 
sion spectroscopy (ARPES) experiments on BaFc2As2 
[3^ show that the band-width is strongly rcnormalized, 
the larger Fermi velocities in electron and hole pockets 
are v e ~ Vh ~ 0.5 eV-A. These values would give c ~ 0.29 
eV-A. Using smaller Fermi velocities would yield half of 
the larger values or c ~ 0.15 eV-A. Within the local mo- 
ment effective J\ a — Ju, — J 2 model, the spin wave velocity 
is given by c = (J\ a + 2 J 2 )a/2 eV-A. From our measured 
Jia + 2 J 2 = 100 ± 20 meV, c ~ 0.28 eV-A which is also 
fairly close to the ARPES results. Therefore, our present 
data do not allow an unambiguous distinction between 
localized and itinerant description of the AF order in 
SrFc2As2 in terms of the spin wave velocity. 

In summary, we carried inelastic neutron scattering ex- 
periments to study low energy spin wave excitations in 
SrFe2As2- The low-temperature spectrum consists of a 
Bragg peak, a spin gap, and sharp spin wave excitations 
at higher energies. Using a simple Heisenberg Hamilto- 
nian, we find J la + 2J 2 = 100 ± 20 meV, J z = 5 ± 1 meV, 
and J s = 0.015 ± 0.005 meV. On warming across T/v, 
there is weak critical scattering and spin-spin correlations 



in the AF wave vector region explored in the paramag- 
netic state, different from the paramagnetic SDW exci- 
tations in Cr. These results are consistent with the AF 
phase transition in SrFe2As2 being first order in nature. 
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